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We present results obtained in QCD with two flavors of non-perturbatively improved Wilson fermions at finite 
temperature on 16^ x 8 and 24"^ x 10 lattices. We determine the transition temperature in the range of quark masses 
0.6 < rriTr/mp < 0.8 at lattice spacing a«0.1 fm and extrapolate the transition temperature to the continuum and 
to the chiral limits. 



1. INTRODUCTION 

In order to obtain predictions for the real world 
from lattice QCD, we have to extrapolate the lat- 
tice data to the continuum and to the chiral lim- 
its. Recently the Bielefeld group P and the CP- 
PACS collaboration 12 using different fermion 
actions obtained consistent values for the criti- 
cal temperature Tc in the chiral limit, albeit on 
rather coarse lattices at A'^f = 4 and 6. Edwards 
and Heller [3j determined Tc for A^f = 4, 6 us- 
ing nonperturbatively improved Wilson fermions. 
We compute Tc on finer lattices with Nt — 8 and 
10 with high statistics. Our results for Nt = 8 
were reported in Ref. 0]. 

2. SIMULATION 

We use fermionic action for non-perturbatively 
improved Wilson fermions: 



various k. The values of k and the corresponding 
number of configurations for 16^ x 8 and 24^ x 
10 lattices can be found in Ref. 4 and Table 1, 
respectively. The number of configurations for 
24^ X 10 lattice is not large enough and results 
for this lattice are preliminary. We use results 
obtained at T=0 to fix the scale. The contour 
plot of lines of constant ro/a and niT^/mp is 
shown in Fig. ^ 
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where Sp is the original Wilson action, Csw was 
calculated in 0|. 

Configurations are generated on 16"^ x 8 {(3 — 
5.2 and 5.25) and 24^ x 10 (/3 = 5.2) lattices at 

*Talk given by Y. N. at Lattice'04. 



Table l:Simulation statistics on 24^^ x 10. 



3. CRITICAL TEMPERATURE 

We use the Polyakov loop susceptibility to de- 
termine the transition point. The critical value 
of K turns out to be Kt=0.1354. Using ro=0.5 fm 
and interpolating ro/a to the critical point, we 
obtain for the crtical temperature: 

Tc = 196(4)MeV, m^/mp = 0.64(3) 
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Figure 1. The lines of constant r^/a and niT^/mp 
at r = 0. Crosses correspond to parameters used 
in simulations of 24^ x 10 lattice. 
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Figure 2. Polyakov loop (top) and its susceptibil- 
ity (bottom) at /? = 5.2 on 24^ x 10. 



4. CONTINUUM AND CHIRAL LIMITS 

At small enough lattice spacing and quark mass 
one can extrapolate the critical temperature Tc to 
the continuum and the chiral limits using formula: 

Tero = T^^^--\^ + Ca{-f + C,{- - , 

ro K Kc 



where y™"!'""* corresponds to the extrapolated 
value of the critical temperature and /3 and 6 are 
critical indices. 

We make an attempt to fit four values for TcTq 
(see Table 3), obtained at rather large quark 
masses, to estimate the parameters in this extrap- 
olation expression. 
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Table 3: Available data for TcTq. 
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Table 4:Fitting results. 

We extrapolate the value of the critical tem- 
perature using different values of 0.54 and 1 as 
1/(36. If the transition in two-flavor QCD is sec- 
ond order, the transition is expected to belong to 
the universality class of the 3D 0(4) spin model 
with 1//3i5r:!0.54. If the transition is first order, 
then 1/(36^1. Table 4 and Figs. H present fit- 
ting results. We get the critical temperature in 
the continuum and in the chiral limits. 



In the case of l//3(5=0.54: 



(1) 



This value agrees with values obtained in Refs.^ 



In the case of 1//3S—1: 
ym„a^o ^ 201(12) MeV. 



(2) 



Although some lattice studies jll2j indicate sec- 
ond order chiral transition in two- flavor QCD, 
there are also results |7| supporting first order 
transition. Results of our fits do not allow to dis- 
criminate between first and second order transi- 
tions because of rather large errors in TcVq values. 
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5. CONCLUSIONS 

We determined the critical temperature in full 
QCD on X 10 lattice at /? = 5.2 with Nf = 2 
clover fermions using Polyakov loop susceptibil- 
ity. Our results are in agreement with the results 
of other groups, as it is shown in Fig. |31 We ex- 
trapolate the critical temperature Tc to the con- 
tinuum and the chiral limits both in cases of first 
and second order transition. The extrapolation 
results are given by and 10). We are contin- 
uing simulations on 24^ x 10 lattice to get better 
precision of value on this lattice. 
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Figure 3. Tcl\fo as a function of m-nlmp. The 
circle and squares show results of jS| and P re- 
spectively. The diamond and triangles corre- 
spond to our data. 
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Figure 4. TcTq as a function of lattice spacing and 
quark mass if the second order phase transition 
is realized. Best fit is shown by the plane. 
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